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ABSTRACT 

It  has  been  predicted  that  point  source  radiation  which 
is  at  near-grazing  incidence  to  a slightly  rough  surface 
generates  a boundary  wave  in  the  fluid  as  well  as  the 
better  known  volume  wave  (I.  Tolstoy,  J.  Acoust.  Soc.  Am., 
63,  s6o  (1978))*  Anechoic  chamber  experiments,  using  a 
surface  constructed  of  rigid  hemispheres  on  a rigid  plate, 
have  been  performed  to  compare  with  theory.  The  experiments 
confirm  that  the  relative  boundary  wave  amplitude  to  volume 
wave  amplitude  increases  as  the  3/2  power  of  the  sound  fre- 
quency when  the  source  and  receiver  are  on  the  boundary. 
However  the  predicted  r^  growth  with  range  reaches  a limit 
at  kr** 100 . For  the  receiver  above  the  rough  surface  the 
experimental  values  of  the  scattered  sound  show  some 
agreement  with  theoretical  predictions. 


4 

J 


1 


TABLE  OF  CONTENTS 


ACKNOWLEDGEMENTS  10 

I.  INTRODUCTION  11 

II.  THEORY  12 

III.  RESEARCH  FACILITIES  15 

A.  OCEAN  ACOUSTICS  LABORATORY  AND  ANECHOIC 

CHAMBER  FACILITIES  15 

B.  EQUIPMENT  LIST  17 

IV.  EXPERIMENTAL  DESIGN  AND  PROCEDURES  19 

A.  TOLSTOY  ROUGH  SURFACE  MODEL  19 

B.  SOUND  SOURCE  AND  RECEIVER  SELECTION  19 

C.  DATA  COLLECTION  PROCEDURES  26 

D.  COMPUTER  PROGRAM  PROCESSES  28 

V.  EXPERIMENTAL  RESULTS  AND  ANALYSIS  36 

A.  SOURCE  AND  RECEIVER  ON  SURFACE  36 

B.  SOURCE  ON  SURFACE,  RECEIVER  ABOVE  SURFACE  - 47 

VI.  CONCLUSIONS  53 

APPENDIX  A 54 

BIBLIOGRAPHY  109 

INITIAL  DISTRIBUTION  LIST  110 


5 


LIST  OF  FIGURES 


FIGURE 

1 

OCEAN  ACOUSTICS  LABORATORY  COMPUTER 

SYSTEM  

16 

FIGURE 

2 

ROUGH  SURFACE  MODEL  DRAWING  

20 

FIGURE 

3 

MAXIMUM  PACKING  DENSITY  OF  THE  HEMISPHERICAL 
BOSSES  - 

21 

FIGURE 

4 

RUBBER  CEMENT  FILL  TO  EQUATORIAL  PLANE 

OF  BOSSES  

21 

FIGURE 

5 

SMOOTH  SURFACE  CONTROL  PLATE  WITH  B & K 

4134  MICROPHONE  POSITIONED  AT  TWENTY 
CENTIMETERS  

22 

FIGURE 

6 

ROUGH  SURFACE  PLATE  WITH  B & K 4134  MIC- 
ROPHONE POSITIONED  AT  TWENTY  CENTIMETERS  - 

23 

FIGURE 

7 

EXPERIMENTAL  PLATES  IN  THE  ANECHOIC 

CHAMBER  

24 

FIGURE 

8 

OCEAN  ACOUSTICS  LABORATORY  RECEIVING  AND 
PROCESSING  SYSTEM  

29 

FIGURE 

9 

SCHEMATIC  OF  THE  OCEAN  ACOUSTICS 

LABORATORY  EQUIPMENT  

30 

FIGURE 

10 

TRANSMITTING  AND  RECEIVING  EQUIPMENT  IN 

THE  ANECHOIC  CHAMBER  

31 

FIGURE 

11 

SCHEMATIC  DRAWING  OF  THE  ANECHOIC  CHAMBER 
EQUIPMENT  

32 

FIGURE 

12 

ONE  CYCLE  OF  THE  2.5  KILOHERTZ  SAWTOOTH 
PULSE.  RECEIVED  PULSE.  THE  SAMPLING 

PULSE.  SAMPLING  APERTURE  PULSE  

33 

FIGURE 

13 

CROSS-SECTIONAL  DRAWING  OF  THE  ROUGH 

SURFACE  MODEL  

33 

FIGURE 

14 

COMPOSITE  PLOT  (RATIO  OF  BWA  TO  VWA  VS. 
FREQUENCY)  10  CM  

41 

FIGURE 

15 

COMPOSITE  PLOT  (RATIO  OF  BWA  TO  VWA  VS. 
FREQUENCY  ) 20  CM  

42 

FIGURE 

16 

COMPOSITE  PLOT  (RATIO  OF  BWA  TO  VWA  VS. 
FREQUENCY)  30  CM  

43 

6 


I 


FIGURE  17 

FIGURE  18 

FIGURE  19 
FIGURE  20 
FIGURE  21 
FIGURE  22 
FIGURE  23 

FIGURE  24 
FIGURE  25 
FIGURE  26 
FIGURE  27 
FIGURE  28 
FIGURE  29 
FIGURE  30 
FIGURE  31 
FIGURE  32 
FIGURE  33 
FIGURE  34 
FIGURE  35 
FIGURE  36 
FIGURE  37 

FIGURE  38 


COMPOSITE  PLOT  (RATIO  OF  BWA  TO  VWA  VS. 
FREQUENCY)  40  CM  

COMPOSITE  PLOT  (RATIO  OF  BWA  TO  VWA  VS. 
FREQUENCY)  50  CM  

bwa/vwa  vs.  range  COMPOSITE  PLOT  

RATIO  OF  Ps  TO  Pg-  VS.  HEIGHT  

RATIO  OF  Ps  TO  Pj.  VS.  HEIGHT  

RATIO  OF  Ps  TO  Pj.  VS.  HEIGHT  

RATIO  OF  Ps  TO  Pg  VS.  HEIGHT  

FIGURES  IN  APPENDIX 

BIOT-TOLSTOY  THEORETICAL  VALUES  

FAST  FOURIER  TRANSFORM  DATA  FORMAT  


COMPUTER  PROCESSED  DATA  FORMAT 


RATIO 

OF 

BWA 

TO 

VWA 

VS. 

FREQUENCY 

(10 

CM) 

RATIO 

OF 

BWA 

TO 

VWA 

VS. 

FREQUENCY 

(20 

CM) 

RATIO 

OF 

BWA 

TO 

VWA 

VS. 

FREQUENCY 

(30 

CM) 

RATIO 

OF 

BWA 

TO 

VWA 

VS. 

FREQUENCY 

(40 

CM) 

RATIO 

OF 

BWA 

TO 

VWA 

VS. 

FREQUENCY 

(10 

CM) 

RATIO 

OF 

BWA 

TO 

VWA 

VS. 

FREQUENCY 

(20 

CM) 

RATIO 

OF 

BWA 

TO 

VWA 

VS. 

FREQUENCY 

(30 

CM) 

RATIO 

OF 

BWA 

TO 

VWA 

VS. 

RANGE  (10  KHZ) 

— 

RATIO 

OF 

BWA 

TO 

VWA 

VS. 

RANGE  (20JCHZ) 

— 

RATIO 

OF 

BWA 

TO 

VWA 

VS. 

RANGE  (30  KHZ) 

— 

FFT  DATA  RELATIVE  AMPLITUDE  VS.  FREQUENCY 
FOR  SMOOTH  PLATE  

FFT  DATA  RELATIVE  AMPLITUDE  VS.  FREQUENCY 
FOR  ROUGH  PLATE  


FIGURE  39  FFT  DATA  RELATIVE  AMPLITUDE  VS.  FREQUENCY 

FOR  HEMISPHERICAL  BOSS  REPACKS  70 

FIGURE  40  RATIO  OF  BWA  TO  VWA  VS.  FREQUENCY  (10  CM)-  71 

FIGURE  41  RATIO  OF  BWA  TO  VWA  VS.  FREQUENCY  (20  CM)-  72 

FIGURE  42  RATIO  OF  BWA  TO  VWA  VS.  FREQUENCY  (30  CM)-  73 

FIGURE  43  RATIO  OF  BWA  TO  VWA  VS.  FREQUENCY  (40  CM)-  74 

FIGURE  44  RATIO  OF  BWA  TO  VWA  VS.  FREQUENCY  (50  CM)-  75 

FIGURE  45  RATIO  OF  BWA  TO  VWA  VS.  RANGE  (10  KHZ)  — 76 

FIGURE  46  RATIO  OF  BWA  TO  VWA  VS.  RANGE  (20  KHZ)  — 77 

FIGURE  47  RATIO  OF  BWA  TO  VWA  VS.  RANGE  (30  KHZ)  — 78 

FIGURE  48  RATIO  OF  BWA  TO  VWA  VS.  FREQUENCY  (10  CM)-  79 

FIGURE  49  RATIO  OF  BWA  TO  VWA  VS.  FREQUENCY  (20  CM)-  80 

FIGURE  50  RATIO  OF  BWA  TO  VWA  VS.  FREQUENCY  (30  CM)-  81 

FIGURE  51  RATIO  OF  BWA  TO  VWA  VS.  FREQUENCY  (40  CM)-  82 

FIGURE  52  RATIO  OF  BWA  TO  VWA  VS.  FREQUENCY  (50  CM)-  83 

FIGURE  53  RATIO  OF  BWA  TO  VWA  VS.  RANGE  (10  KHZ)  — 84 

FIGURE  54  RATIO  OF  BWA  TO  VWA  VS.  RANGE  (12.5  KHZ)—  85 

FIGURE  55  RATIO  OF  BWA  TO  VWA  VS.  RANGE  (15  KHZ)  — 86 

FIGURE  56  RATIO  OF  BWA  TO  VWA  VS.  RANGE  (17-5  KHZ)—  87 

FIGURE  57  RATIO  OF  BWA  TO  VWA  VS.  RANGE  (20  KHZ) 88 

FIGURE  58  RATIO  OF  BWA  TO  VWA  VS.  RANGE  (22.5  KHZ)—  89 

FIGURE  59  RATIO  OF  BWA  TO  VWA  VS.  RANGE  (25  KHZ) 80 

FIGURE  60  RATIO  OF  BWA  TO  VWA  VS.  RANGE  (27-5  KHZ)—  91 

FIGURE  61  RATIO  OF  BWA  TO  VWA  VS.  FREQUENCY  (10  CM)-  92 

FIGURE  62  RATIO  OF  BWA  TO  VWA  VS.  FREQUENCY  (20  CM)-  93 

FIGURE  63  RATIO  OF  BWA  TO  VWA  VS.  FREQUENCY  (30  CM)-  94 


8 


FIGURE  64  RATIO  OF  BWA  TO  VWA  VS. 
FIGURE  65  RATIO  OF  BWA  TO  VWA  VS. 
FIGURE  66  RATIO  OF  BWA  TO  VWA  VS. 
FIGURE  67  RATIO  OF  BWA  TO  VWA  VS. 
FIGURE  68  RATIO  OF  BWA  TO  VWA  VS. 
FIGURE  69  RATIO  OF  BWA  TO  VWA  VS. 
FIGURE  70  RATIO  OF  BWA  TO  VWA  VS. 
FIGURE  71  RATIO  OF  BWA  TO  VWA  VS. 
FIGURE  72  RATIO  OF  BWA  TO  VWA  VS. 
FIGURE  73  RATIO  OF  BWA  TO  VWA  VS. 
FIGURE  74  RATIO  OF  BWA  TO  VWA  VS. 


FREQUENCY  (10  CM)-  96 
FREQUENCY  (20  CM)-  97 
FREQUENCY  (30  CM)-  98 
FREQUENCY  (40  CM)-  99 
FREQUENCY  (50  CM)-  100 
RANGE  (10  KHZ)  — 101 

RANGE  (12.5  KHZ)—  102 

RANGE  (15  KHZ) 103 

RANGE  (17.5  KHZ)—  104 
RANGE  (20  KHZ) 105 


FIGURE  75  RATIO  OF  BWA  TO  VWA  VS.  RANGE  (22.5  KHZ)—  106 

FIGURE  76  RATIO  OF  BWA  TO  VWA  VS.  RANGE  (25  KHZ) 107 

FIGURE  77  RATIO  OF  BWA  TO  VWA  VS.  RANGE  (27-5  KHZ)—  108 


9 


ACKNOWLEDGEMENTS 


The  writer  wishes  to  express  his  appreciation  to 
Professor  Herman  Medwin  of  the  Physics  Department,  U.S. 
Naval  Postgraduate  School  and  Professor  Ivan  Tolstoy, 
Visiting  Distinguished  Professor,  U.S.  Naval  Postgraduate 
School  for  their  guidance  and  encouragement  during  the 
preparation  of  this  thesis;  to  Mrs.  Jeanie  Savage  for  her 
computer  programming;  to  Mr.  Tom  Maris  for  his  construc- 
tion of  the  experimental  surfaces  and  the  associated  hard- 
ware; and  especially  to  Lieutenant  Commander  Joe  Bremhorst 
for  his  electronic  wizardry  and  friendship. 


10 


I.  INTRODUCTION 


1 


At  the  95'th  meeting  of  the  Acoustical  Society  of  America 
(May  1978),  I.  Tolstoy  outlined  his  application  of  a theory 
by  M.A.  Biot  describing  sound  scatter  from  slightly  rough 
surfaces  at  near-grazing  incidence.  This  theory,  using 
a boundary  condition  due  to  Biot  (Biot  1968) , is  first 
order  in  the  roughness  parameter,  rather  than  second  order 
as  other  scattering  theories:  it  includes  multiple  scatter 
and  diffraction  which  other  theories  generally  ignore:  it 
does  not  involve  the  Kirchhoff  Assumption  which  is  partic- 
ularly suspect  at  small  grazing  angles. 

For  near-grazing  forward  scatter  from  a point  source, 

« 

the  theory  predicts  that,  in  addition  to  a volume  wave,  co- 
herent  multiple  Rayleigh  scatter  will  generate  a boundary 
wave  in  the  fluid  above  the  slightly  rough  rigid  surface. 

At  sufficient  distances  for  near-grazing  incidence  it  is 
predicted  that  the  boundary  wave  will  be  greater  than  the 
direct  or  volume  wave . 

The  aim  of  this  work,  in  general,  is  to  compare  Tolstoy's 
theory  with  experimental  data. 
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II.  THEORY 


Tolstoy,  using  an  original  formulation  of  coherent 
scatter  proposed  by  Biot  (1968,  1973).  has  extended  the 
theory  to  obtain  solutions  for  the  scattering  of  transient 
spherical  waves  by  rough  planes,  providing  the  roughness 
scale  is  small  compared  to  the  pulse  width.  Biot  had  demon- 
strated that  if  the  spacing  (h)  between  hemispherical  boss 
centers  was  small  compared  to  the  wavelength,  then  these 
bosses  co.'ld  be  replaced  by  a continuous  distribution  of 
monopoles  and  dipoles.  Beginning  with  the  assumption 
&/\<h/\«l , where  d is  a radius  of  a boss  and  h is  the 

center  to  center  spacing  of  the  bosses,  he  proceeded  to 

« 

.derive  the  equations  for  the  scattering  of  pressure  pulses 
radiated  by  transient  point  sources  above  or  on  a rough 
plane  (z=0)  for  an  isotropic  case.  After  extensive  math- 
ematics, beyond  the  scope  of  this  thesis,  Tolstoy  arrived 
at  a temporal  expression  for  the  boundary  wave  pressure. 
Tolstoy  also  carried  his  solution  to  the  frequency  domain 
with  the  use  of  Fourier  transform  techniques.  It  is  here, 
from  Tolstoy's  equation  for  an  impulse  Fourier  transform 
for  the  boundary  wave  solution,  that  this  experimental 
verification  starts. 

The  spectral  amplitude  of  the  boundary  wave  (fjJ)  given 
by  Tolstoy  (1979)  is 


~C* 


i 


l 
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where  ^ is  a scattering  parameter  (proportional  to  volume 
of  scatterers  per  unit  area) , and 


Jii_  = k 

C K‘ 


For 


kr»  1 , 

far  field  substitutions  for  the  Bessel  functions  allow  the 
equation  for  the  total  scattered  field  to  be  written  as: 


Pstfl-£k' 


’<*♦£  2,2ejV  Sm(^-X)* 


(2) 


where  for  close-packed  hemi- 
spherical bosses  of  diameter 
0.2  cm  in  air 
£,  = 8.85xl0~^cm 


-1 


k = w/c  = 2irf  ..cm 

J^xlO4 


r = range  of  source  to  receiver  in  cm 
f = frequency  in  Hertz 
zQ=  height  of  source  in  cm 
z = height  of  receiver  in  cm 
•©■=  tan"1  -z— -Q 


For  the  case  of  the  source  and  receiver  on  the  plate 
(zQ  + z = 0)  this  equation  reduces  to 

-1/  3/ 


Ps<f)  = petf ) - ^ k1  ]^r]  2 - ^ (Z.tt  r)  2 k 


(3) 


which  can  be  readily  solved  for  various  range  and  frequency 
combinations.  The  direct  volume  acoustic  arrival,  for  the 
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impulsive  case  is  $(t  - r/c)  and  lts  spectrum. 

IJ.  = [2rrrj' 


-I 


(4) 


It  is  the  ratio 


Pa/P*  =e(2irr)'^  k3/i 


for  z„  = z = 0 
o 


(5) 


calculated  for  various  ranges  and  frequencies,  that  is 
compared  with  experimental  spectral  amplitudes  as  the  major 
effort  in  this  thesis.  (See  Section  IV  - D for  experimental 
procedures . ) 
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A.  OCEAN  ACOUSTICS  LABORATORY  AND  ANECHOIC  CHAMBER 

Due  to  the  close  proximity  of  the  Ocean  Acoustics 
Laboratory,  where  the  data  collection  and  processing  equip- 
ment was  located,  and  the  Anechoic  Chamber,  where  the  ex- 
perimental apparatus  was  located,  it  was  possible  to 
collect  the  analog  signal  and  digitize  it  in  real  time. 

Data  acquisition  and  processing  were  accomplished  using  a 
digital  computer  system  composed  of  four  primary  components, 
interfaced  to  provide  high  speed  analog  to  digital  conver- 
sion, digital  processing,  data  storage  and  data  printout. 

The  design  was  developed  by  the  Special  Projects  Section 
of  the  Naval  Air  Development  Center  in  conjunction  with 
Pinkerton  Computer  Consultants,  Inc.  of  Warminster, 
Pennsylvania.  The  four  individual  components  are  described 
belowi  (see  Figure  1) 

1 . Interdata  Model  70  Computer 

This  digital  minicomputer  is  FORTRAN  and  BASIC 
programmable  with  a 64  thousand  byte  core  memory.  In  ad- 
dition to  the  core  memory,  data  can  be  stored  on  magnetic 
discs . 

2.  Phoenix  Analog  to  Digital  Converter,  Model  ADC  712 
The  ADC  analog  to  digital  converter  is  a high  speed 

device  capable  of  encoding  +10  volt  input  signals  in  digital 
form  and  providing  an  accuracy  of  0.005  volts  in  20  volts. 
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The  maximum  sampling  rate  is  450,000  samples  per  seconds 
however,  rates  of  256,000  samples  per  second  and  320,000 
samples  per  second  were  used  in  the  experiment  for  optimum 
experimental  frequency  resolution. 

3-  Texas  Instruments  Silent  Electronic  Data  Terminal, 
Model  733 

The  TI  733  consists  of  a keyboard  used  as  a program- 
ming input/output  control  device,  a printer,  and  a playback/ 
record  section  used  in  conjunction  with  digital  cassettes. 

4.  Orbis  Model  76  Diskette  Drive 

The  Orbis  Model  76  Diskette  Drive  is  a small,  port- 
able, direct  access,  256  kilobyte  floppy  disc  data  storage 
device.  These  magnetic  discs  provide  the  capability  to 
store  large  quantities  of  data  for  later  computer  analysis. 
The  overall  system  facilitates  rapid,  accurate  processing 
of  any  desired  type  of  analog  electrical  signal  and  was  used 
mainly  for  frequency  analysis  using  standard  Fast  Fourier 
Transform  algorithms. 


B . EQUIPMENT  LIST 

A list  of  standard  equipment  referred  to  throughout 

the  text  is  given  belowt 

Abbreviation  Full  Description 

Scope  Tektronix  Model  545B  4-channel 

Oscilloscope 

HP  Counter  5223L  Hewlett-Packard  Model  5223L 

Electronic  Counter 

1 


RS648  Timing 
Simulator 

Wavetek  114 

Wavetek  114 

Frequency  Counter 


RMS  Voltmeter 


Interface  Technology  Model  RS648 
Timing  Simulator 

Wavetek  Model  114  Function  Generator 

Wavetek  Model  144  Function  Generator 

Donner  Model  8050  Frequency/Period 
Meter 

Fluke  8920A  RIVE  Voltmeter 


Digital  Multimeter  Fluke  8000A  Digital  Multimeter 


Pulse  Transformer 


DC  Power  Supply 
KH  3550  Filter 


PAR  133  Preamp 


MIC  Power  Supply 


Source 

Receiver 

O-scope 


MIC  Preamp 


North  Hills  Wide-Band  Pulse 
Transformer 

Lambda  Regulated  Power  Supply 

Krohn-Hite  Model  3550  Frequency 
Filter 

Princton  Applied  Research  Model  113 
Preamp 

B&K  Model  2804  Microphone  Power 
Supply 

B&K  Model  4134  Microphone 

B&K  Model  4133  Microphone 

Hewlett-Packard  Model  140A  Dual 
Trace  Oscilloscope 

B&K  Model  2619 


IV.  EXPERIMENTAL  DESIGN  AND  PROCEDURES 


A.  TOLSTOY  ROUGH  SURFACE  MODEL 

The  experiment  required  that  sound  at  near-grazing  in- 
cidence be  sent  and  received  over  a slightly  rough  rigid 
surface.  More  specifically,  the  surface  was  to  be  composed 
of  hemispherical  bosses  with  maximum  hexagonal  packing 
density.  To  achieve  these  parameters  two  rectangular  pieces 
of  aluminum  plate  were  cut  with  the  dimension  of  100  cen- 

| 

timeters  by  82  centimeters.  One-half  inch  holes  were  drill- 
ed and  tapped  along  the  center  line  of  each  plate  at  10  cen- 
timeter intervals.  Threaded  plugs  and  two  non-conducting 
threaded  sleeves  were  machined  to  hold  the  source  and 
receiver  at  various  distance  combinations  on  the  plates. 

Next,  edges  were  fitted  on  one  plate  and  a layer  of  #9  chill- 
ed lead  shot,  2 millimeters  in  diameter  each,  was  poured  on 
the  surface.  Using  shaker  table  methods  on  the  slightly  el- 
evated plate,  maximum  packing  was  achieved.  (See  Figure  3) 
Thinned  rubber  cement  was  then  poured  over  the  shot  to  in- 
sure adherence  to  the  plate  and  to  fill  the  spaces  up  to  the 
equator  of  the  spheres.  (See  Figure  4)  The  smooth  and  the 
rough  surface  plates  were  then  suspended  by  wires  in  the 
anechoic  chamber.  (See  Figures  5*  6,  and  7) 

B.  SOUND  SOURCE  AND  RECEIVER  SELECTION 

Several  sources  were  investigated  including  an  Ionovac 
exponential  horn,  a B&K  Model  4138  (1/8  inch)  condenser 


♦Source  located  in  this  position  for  all  data  presented 
unless  otherwise  specified. 

Figure  2.  Rough  Surface  Model  Drawing 


20 


Figure  3-  Maximum  Packing  Density  Of  The  Hemispherical 
Bosses  (Magnification  15x) 


Figure  4.  Rubber  Cement  Fill  to  Equatorial  Plane  Of 
Bosses  (Magnification  33x) 
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Smooth  Surface  Control  Plate  with  B + K 413 
Microphone  Positioned  at  Twenty  Centimeters 
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Lgure  . Experimental  Plates  Tn  The  Anechoic  Chamber 


microphone,  a B&K  Model  4136  (1/4  inch)  condenser  microphone, 
a B&K  Model  4133  and  a Model  4134  (£  inch)  condenser  micro- 
phone and  a B&K  Model  4145  (1  inch)  condenser  microphone. 

With  these  piston-type  sources  embedded  flush  in  the 
plate  and  by  pulsing  the  sent  signal,  it  was  possible  to  pro- 
duce acoustic  pressures  almost  identical  with  those  produced 
by  a hemispherical  radiating  simple  (point)  source  of  equal 
strength,  except  for  the  presence  of  a directivity  term. 

The  directivity  was  investigated  for  various  frequencies  as 
a function  of  source  radii  (l/l6,  l/8,  l/4,  l/2  inch).  The 
1/8  inch  diameter  B&K  microphone  most  closely  approximated 
a point  source,  however  the  maximum  acoustic  pressure  gener- 
ated was  lower  than  for  larger  sources.  The  B&K  Model  4134 
(1/2  inch)  microphone  was  selected  as  the  best  compromise 
for  large  output  with  nearly  omnidirectional  radiation  pat- 
tern. This  microphone  could  be  driven  with  a maximum  of 
250  volts  DC  and  35  volts  (RMS)  AC  and  approximated  adequate- 
ly a point  source  for  the  wave-lengths  to  be  used. 

Receiver  selection,  likewise,  was  a compromise  based 
on  the  sensitivity  of  the  receiver  and  its  frequency  re- 
sponse range.  The  B&K  Model  4133,  with  an  8 centimeter  probe 
tube  of  outer  diameter  3.2  millimeters  was  selected  as  the 
receiver.  With  a frequency  range  from  3.9  Hertz  to  40  kilo- 
hertz and  a noise  floor  at  -38  dB  re  1 Pa,  this  receiver 
provided  the  very  best  sensitivity  and  lowest  noise  available 
A compatible  B&K  Model  2804  DC  microphone  power  supply  was 
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used  in  conjunction  with  a B&K  Model  2619  preamplifier  for 
this  l/2  inch  condenser  microphone.  This  power  supply  was 
electrically  quieter  and,  of  prime  importance,  did  not  drive 
the  microphone  heating  element,  thus  reducing  thermally 
induced  signal  fluctuations. 

C.  DATA  COLLECTION  PROCEDURES 

The  experiment,  in  general,  involved  the  positioning  of 
a source  and  receiver  alternately  on  the  smooth  and  the 
rough  surface  plates  and  recording  the  received  signal  for 
subsequent  computer  evaluation.  As  can  he  seen  from  figures 
2 and  5>  "the  source  and  receiver  could  he  located  in  any  of 
nine  sites  on  either  plate.  These  positions  as  shown  in 
figure  2,  will  he  referred  to  in  the  remainder  of  the  text 
hy  their  xn  designator.  These  source  and  receiver  points 
were  located  at  ten  centimeter  intervals,  center  to  center. 

With  the  source  and  receiver  positioned,  a 2.5  kHz  pulsed 
sawtooth  was  sent  hy  the  Wavetek  Model  144  through  a wide- 
band pulse  transformer  with  a 4 to  1 step-up  voltage.  An 
example  of  the  sent  signal  can  be  seen  in  the  top  trace  of 
figure  12 . The  slight  distortion  in  the  sawtooth  was  a 
result  of  the  transformer  response  to  this  relatively  low 
frequency  signal.  However  a flat  response  was  not  essential 
to  the  experiment.  The  output  AC  voltage  was  maintained  at 
a constant  35  volts  RMS.  A 250  volt  DC  polarizing  voltage 
was  then  superimposed  for  a total  of  300  volts  peak  applied 


to  the  B&K  Model  4134  source.  These  voltages  and  the  source 
frequency  were  continually  monitored  for  stability  using  a 
Donner  Frequency  Counter,  a Fluke  True  RMS  Voltmeter,  a 
Fluke  Digital  Multimeter  and  a Hewlett-Packard  Dual  Trace 
Oscilloscope  (see  Figures  10  and  11). 

The  receiver  signal  (see  second  trace  of  Figure  12)  from 
the  B&K  Model  4133  microphone,  with  its  associated  2619  Pre- 
amplifier and  2804  Power  Supply,  was  then  fed  through  a 
PAR  133  PreAmp  with  a gain  of  10.  It  then  passed  through 
two  cascaded  Krohn-Hite  Model  3550  Filters,  high  passing 
above  1 kHz  and  low  passing  below  100  kHz.  The  filtered 
signal  was  then  amplified  by  another  PAR  113  PreAmp  with  a 
gain  of  200.  • 

The  received  analog  signal  was  then  sent  into  the  ad- 
joining Ocean  Acoustics  Laboratory  to  a Phoenix  Analog  to 
Digital  Converter  and  digitized.  Generally  1000  blocks  of 
this  pulsed  signal  were  then  averaged  before  being  Fast 
Fourier  Transformed  by  standard  FFT  algorithms.  The  results 
were  then  stored  on  magnetic  discs. 

In  the  Ocean  Acoustics  Laboratory,  an  Interface  Tech- 
nology Model  RS648  Timing  Simulator  was  used  as  the  master 
timer  for  the  sawtooth  generator,  the  sampling  frequency  gen- 
erator and  the  actual  sampling  duration  of  the  Interdata  70 
Computer.  Using  available  channels  on  the  timing  simulator 
it  was  possible  to  trigger  the  sent  signal,  the  sampling 
frequency  and  the  proper  sampling  aperture  from  the  master 


clock  within  100  nanosecond  accuracy  to  sample  any  analog 
signal  of  interest.  (See  Figures  8 and  9-) 


With  a pulse  duration  of  0.4  milliseconds,  there  was  a 
frequency  resolution  of  2.5  kilohertz.  The  sampling  fre- 
quency, was  in  general,  320  kilohertz  which  provided  128 
samples  in  each  block.  One  thousand  blocks,  each  containing 
128  digitized  samples  were  time  averaged  and  these  values 
printed  out  and/or  stored  in  preparation  for  the  FFT. 

D.  COMPUTER  PROGRAM  PROCESSES 

Theoretical  values  were  calculated  from  equation  (5)  for 
ranges  from  10  centimeters  to  50  centimeters  and  for  fre- 
quencies from  10  kilohertz  to  30  kilohertz.  (See  Figure  24, 
Appendix  A.)  Plotting  these  values,  as  a function  of  fre- 
quency or  range  for  a specific  frequency,  provided  the 
Tolstoy  theoretical  curves. 

The  processing  of  the  experimental  data  for  comparison 
with  theory  was  accomplished  in  the  following  manner.  Sig- 
nals from  the  smooth  plate  and  rough  surface  plate  for  equiv- 
alent distances  were  taken,  Fourier  transformed  and  stored  on 
magnetic  discs.  These  same  data  were  also  printed  out  at 
the  terminal  for  real-time  monitoring  for  consistency  and 
stability.  The  data  were  printed  out  for  frequencies  from 
2.5  kilohertz  to  60  kilohertz  in  individual  columns  titled 
amplitude,  phase  angle,  A (real  part)  and  jB  (imaginary  part) 
of  the  complex  form.  (See  Figures  25  and  26,  Appendix  A for 
an  example  of  the  data  format) . 
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SOURCE  RECEIVER 

Figure  13-  Cross-Sectional  Drawing  Of  The  Rough  Surface 
Model . 


With  the  transformed  data  stored  for  various  distances, 
a computer  program  was  written  to  form  the  required  normal- 
ized ratio  of  the  boundary  wave  amplitude  to  the  volume 
wave  amplitude,  BWA/VWA.  The  program  first  subtracted  the 
complex  spectral  amplitudes  for  the  smooth  plate  (VWA) 
from  the  data  for  the  rough  surface  plate  ; for  the  source 
and  receiver  on  the  surface  the  difference  was  the  complex 
spectral  amplitude  of  the  boundary  wave.  The  real  and  the 
imaginary  components  of  these  complex  boundary  wave  values 
were  then  squared  and  the  square  root  taken  with  the  result 
being  the  boundary  wave  amplitude,  BWA.  The  boundary  wave 
amplitude  for  specific  frequencies  was  then  divided  by  the 
volume  wave  amplitude  (smooth  surface  amplitude)  far  cor- 
responding frequencies  and  their  ratio,  BWA/VWA,  plotted 
versus  frequency  and  range.  These  data  plots  were  then  com- 
pared to  the  theoretical  curves  calculated  from  eg.  5»  (See 
Section  V for  results  and  analysis.) 

For  the  vertical  pressure  evaluations  (Pg/P^.),  the  data 
were  calculated  using  the  same  program  as  for  the  preceding 
process.  However,  equation  (2)  was  now  evaluated  using  in- 
cremental vertical  values  (z)  for  a specific  horizontal  dis- 
tance. Data  were  then  collected  by  positioning  the  receiver 
probe  at  the  equivalent  required  vertical  positions.  The 
experimental  P was  determined  by  subtracting  the  smooth 
surface  complex  spectral  amplitude  from  the  rough  surface 
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complex  spectral  amplitude  and  dividing  this  amplitude  by 
the  smooth  surface  amplitude. 

The  theoretical  curve  of  Pg/P$  was  determined  by 
dividing  equation  (2)  by  the  spectral  amplitude  for  an  ini' 
pulse  in  half  space  (eq.  4). 
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V.  EXPERIMENTAL  RESULTS  AND  ANALYSIS 


A.  SOURCE  AND  RECEIVER  ON  SURFACE 

The  first  rough  surface  model  (Plate  1)  was  constructed 
with  a source/receiver  hold  drilled  in  the  center  of  the 
plate  and  three  equally  spaced  holes  on  either  side,  as 
described  in  Section  IV,  A.  The  source  was  positioned  in 
hole  x^  (20  cm.  from  the  end,  see  Figure  2)  and  the  re- 
ceiver was  then  moved  form  Xgto  x^  and  finally  to  x^.  For 
these  initial  runs  two  cycles  of  the  2.5  kilohertz  sawtooth 
were  sent  and  a sampling  rate  of  256  kilohertz  was  used. 

The  B&K  4133  receiver  was  positioned  flush  (z=0)  in  the 
plates  without  the  microphone  grid  installed.  A run  was 
made  at  10  centimeters  and  every  fourth  point  was  printed 
out  giving  the  desired  2500  Hertz  multiples  for  evaluation. 
This  data  was  compared  favorably  with  previous  recorded  data 
and  showed  excellent  consistency  up  to  approximately 
25  kilohertz. 

In  hopes  of  improving  higher  frequency  stability  (~25 
kilohertz  and  above),  a B&K  O.32  centimeter  outer  diameter 
probe  tube  was  attached  to  the  receiver.  All  subsequent  runs 
were  made  with  tube  installed,  primarily  to  protect  the  ex- 
posed diaphragm,  since  very  little  enhancement  was  observed 
at  the  frequencies  of  interest. 

These  data  showed,  in  general,  excellent  agreement  with 
theory.  (See  Figures  27  through  30»  Appendix  A). 
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Next,  two  runs  were  made  reversing  the  B&K  4133  and 
4134  microphones.  The  data  collection  procedures  were  mod- 
ified by  taking  a data  set  at  a given  distance  on  a given 
plate  and  then  after  several  minutes,  another  data  set  was 
taken  before  the  source  and  receiver  were  changed  to  the 
next  position.  The  first  smooth  data  set  was  then  processed 
with  the  first  rough  data  set  and  then  with  the  second 
rough  data  set.  This  procedure  was  also  used  for  the  second 
smooth  data  set.  The  results  were  then  plotted.  (See 
Figures  31  through  33,  Appendix  A).  The  excellent  agreement 
with  the  theoretical  curves  of  these  data,  would  appear  to 
have  ruled  out  any  adverse  temporal  effects  on  the  experiment. 

The  BWA/VWA  of  these  same  data  were  then  plotted  for 
specific  frequencies  (10  kHz,  20  kHz,  and  30  kHz)  as  a 
function  of  range.  (See  Figures  34  through  36,  Appendix  A.) 
These  preliminary  data  also  showed  excellent  agreement  with 
theory  to  their  limit  of  30  centimeters.  However,  further 
evaluation  was  required  at  greater  source  to  receiver 
ranges . 

Before  continuing  further,  it  was  decided  to  try  and 
determine  if  any  significant  spatially  induced  data  fluc- 
tuations were  generated  by  the  frequent  positioning  and 
subsequent  adjustments  of  the  source  and  receiver  in  the 
plates,  To  check  this,  the  source  was  removed  and  replaced. 
Data  were  collected  and  compared  for  each  individual  run. 

(See  Figures  37  and  38,  Appendix  A for  multiple  data 
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consistency.)  From  this  test  it  was  concluded  that  the 
positioning  of  the  source  and  receiver  at  z=0  could  be 
accomplished  to  within  approximately  a millimeter,  and  thus 
that  the  positioning  could  be  considered  equivalent  for 
all  runs. 

Due  to  the  fact  that  the  source/receiver  positioning 
holes  were  one-half  inch  in  diameter,  it  was  necessary  to 
fill  in  with  dry  bosses  each  time  the  receiver  was  reposi- 
tioned in  the  rough  surface  model.  To  investigate  the  sen- 
sitivity of  this  parameter  on  the  experiment,  five  consecutive 
packings  were  accomplished  without  any  significant  fluctua- 
tions in  the  transformed  data.  Thus  it  was  concluded  that 
the  repacking  required  for  each  positioning  of  the  receiver 
could  also  be  considered  constant.  (See  Figure  39,  Appendix 
A for  data. ) 

Source/receiver  positioning  holes  xQ  and  Xg  were  drilled 
in  both  plates  to  extend  measurement  capabilities  to  maximum 
range  available.  However  aluminum  shavings  adhered  to  the 
surface  of  the  rough  plate  and  a new  surface  had  to  be  con- 
structed (Plate  2) . It  was  also  decided  to  try  a dilute 
polyurethane  as  the  adhesive/filler  vice  the  rubber  cement. 

Two  data  runs  were  made  with  the  source  positioned  at  xQ 
and  the  receiver  progressively  positioned  at  x^  through  x^. 

The  plotted  data  showed  considerable  variation  when  compared 
to  the  previous  results  for  Plate  1.  (See  Figures  40  through 
47,  Appendix  A.) 
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Further  investigation  revealed  that  the  wavelength  of 
the  two  cycles  of  the  sawtooth,  coupled  with  the  minimal 
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20  centimeters  travel  between  the  source  and  edge  of  the 
plate  was  inducing  coherent/incoherent  signal  interference. 

In  addition  microscopic  inspection  of  this  surface  revealed 
that  the  bosses,  in  fact,  were  covered  to  a greater  depth 
than  their  equatorial  planes  and  that  there  were  air  spaces 
under  the  equatorial  plane.  However,  these  data  still 
showed  fair  agreement  with  theory. 

A third  rough  surface  (Plate  3)  was  constructed  with 
rubber  cement  as  the  adhesive/filler . One  cycle  of  the  2.5 
kilohertz  sawtooth  signal  was  sent  with  the  source  position- 
ed at  to  Insure  that  the  previously  stated  problems 
would  be  avoided  completely.  These  data  once  again  showed 
excellent  agreement  with  the  Tolstoy  theoretical  curves  up 
to  the  27-5  to  30  kilohertz  region. 

These  data  were  also  plotted  as  a function  of  range  for 
specific  frequencies  and  compared  with  previous  data  runs. 
They  also  showed  excellent  agreement  with  theory  at  fre- 
quencies up  to  about  15  kilohertz  for  all  tested  ranges. 
However,  as  the  frequency  increased  these  data  diverged 
with  range  relative  to  the  theoretical  curves.  (See  Figures 
48  through  60,  Appendix  A.)  It  was  felt  that  attenuation 
effects  were  the  primary  reason  for  this  divergence,  as  these 
effects  are  not  considered  in  the  theory. 


Since  this  was  the  third  separate  experimental  rough 
surface  and  since  the  results  were  consistent  with  theory 
to  approximately  30  kilohertz,  it  was  decided  to  attempt 
to  intentionally  smooth  the  rough  surface  and  thus  change 
the  roughness  parameter  (4) . Lightweight  lubricating  oil 
was  poured  over  the  bosses  until  only  20-25  percent  of  the 
surface  area  of  individual  shot  remained  exposed.  A data 
run  was  made  and  the  results  plotted  on  Figures  6l  through 
64,  Appendix  A,  as  a function  of  frequency. 

Because  the  roughness  had  been  decreased  (smaller  fe), 
as  expected  these  data  were  displaced  well  below  the  pre- 
vious data  for  the  hemispherical  packing  density.  A set 
of  smooth  data  was  substituted  as  rough  surface  data  and 
processed  to  evaluate  what  happened  experimentally  as  the 
difference  between  the  smooth  (control)  plate  and  the  rough 
surface  was  decreased  (i.e.  roughness  parameter  reduced). 
These  results  showed  the  curve  displaced  well  below  the 
theoretical  curve  for  the  hemispherical  bosses.  In  fact, 
the  ratio  BWA/VWA  differed  from  zero  only  because  of  noise 


in  the  signal  processing.  This  procedure  was  undertaken 
simply  as  a sidelight  to  investigate  the  relative  sensi- 
tivity of  the  experiment  and  theory  to  the  surface  roughness. 
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Composite  Plot  (Ratio  of  BWA  to  VWA  vs.  Frequency) 
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A fourth  plate  was  only  partially  constructed  and  never 
used.  Plate  5 was  the  final  rough  surface  model  constructed. 
A final  data  run  was  made  to  a range  of  50  centimeters  and 
these  data  plotted.  (See  Figures  65  through  77.  Appendix  A.) 
Once  again  they  showed  agreement  with  theory  except  that  at 
increasing  range  the  ratio  BWA/VWA  drops  away  from  the 
theoretical  curve.  This  appears  to  happen  at  krdlOO  and  is 
presumably  due  to  attenuation.  The  composite  Figures  14 
through  19  summarize  the  experimental  work  for  source  and 
receiver  on  surface. 

B.  SOURCE  ON  SURFACE;  RECEIVER  ABOVE  SURFACE 

The  source  and  receiver  were  positioned  twenty  centi- 
meters apart  and  data  sets  were  taken  for  individual  runs 
by  moving  the  receiver  probe  above  the  plates  in  one  centi- 
meter increments  to  five  centimeters.  These  processed  data 
were  then  plotted.  (See  Figures  20  and  21.)  This  entire 
run  showed  much  divergence  when  compared  with  theory,  due 
to  the  fact  that  at  twenty  centimeters  horizontal  distance 
and  for  vertical  distances  above  two  centimeters  the  angle 
of  inclination  was  excessive  to  compare  with  theory. 

The  same  experiment  was  repeated  with  data  sets  being 
collected  at  one  quarter  centimeter  intervals  to  a maximum 
height  of  two  centimeters.  These  data  were  processed  and 
plotted.  The  ratio,  Ps  to  P$  , was  plotted  against  range 
for  specific  frequencies.  The  results  were  then  compared 
to  calculated  theoretical  values.  Fair  agreement  with 


47 


>- 

__ 

cc 

__ 



♦ 

o 

4 

■ 

— — 

UJ 

■ 

* 

N 

X 

M 

M 

X 

X 

X 

* 

m 

o 

lO 

CJ 

(M 

VI.  CONCLUSIONS 


Using  the  experimental  rough  surface  models,  it  has 
been  shown  that  for  kh^l  and  for  ranges  less  than  krslOO, 
when  source  and  receiver  are  on  the  surface  the  near-grazing 
incidence  scattering  is  in  excellent  agreement  with  the 
Tolstoy  predictions  for  the  boundary  wave  behavior.  (See 
Figures  14  through  18.)  As  the  range  is  increased  beyond 
krarlOO  the  experimental  data  diverges  from  theory.  (See 
Figure  19*)  This  divergence  is  due  presumably  to  boundary 
wave  attenuation  processes  which  the  theory  does  not  include. 

For  the  receiver  above  the  surface  the  experimental 
values  of  the  scattered  sound  show  some  agreement  with  the 
theoretical  predictions:  for  example  the  dependence  on 
receiver  height  is  approximately  parallel  to  theory. 
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Figure  25.  Fast  Fourier  Transform  Data  Format. 
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Figure  26.  Computer  Processed  Data  Format 
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Figure  27.  Ratio  of  BWA  to  VWA  vs.  Frequency  (10  cm. 


Figure  28.  Ratio  of  BWA  To  VWA  vs.  Frequency  (20  cm. 
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Figure  3 0.  Ratio  Of  BWA  To  VWA  vs.  Frequency  (40  cm. 


• *3  OCT  78 
PLATE  1 10  CM 


LU  UJ 
CO  CD 
O O 
C£  QC 

a.  a 

QC  DC 

> > £ 
UO  5 

cc  cc  O 
ui 

J 

fO  rO  p 


CM  CM 


VMA  01  VM9  JO  0I1V8 


62 


Figure  31.  Ratio-  Of  BWA  To  VWA  vs.  Frequency  (10  cm. 


• * 3 OCT  78 
PLATE  1 20  CM 


• * 3 OCT  78 
PLATE  1 30  CM 


io  of  BWA  To  VWA  vs.  Frequency  (30  cm. 


• * 3 OCT  78 
PLATE  1 10KHZ 


£ £ 


o O UJ 
> O AJ 

G i-  a- 


™ § < 

1/1  M 

^ 

M ^ CO 

ii  in  m 

N N N 


— 

OJ 

UJ 

UJ 

03 

(D 

O 

cn 

s 

Q. 

Q. 

or. 

CE 

> 

> 

> 

CE 

cc 

C_> 

CE 

o 

UJ 

s 

*T 

rO 

X 

1 — 

• 

4 

1 

VMA  Oi  YM8  dO  OllVd 


vs.  Range  (10  kHz) 


• * 3 OCT  78 
PLATE  1 20  KHZ 


Of  BWA  To  VWA  vs.  Range  (20  kHz) 


| ~ I ■ H ■ " T'  | , . , 


00 

r- 


tn 

r^- 

CM 


r— 

CM 

M" 

< 

< 

fO 

P 

(- 

r~ 

< 

< 

M" 

CD 

Q 

X 

X 

X 

O 

o 

o 

ID 

ID 

3 

O 

O 

O 

CC 

X 

X 

• 

♦ 

♦ 

£ 

1 
+ 
4 

2 

* 


• ♦* 


♦ ♦ 


♦ ♦ • 


M\ 

rO 


O 

rO 


CO 


CM 

. oi 


"oo‘  1 vf7" 


o 

r-- 


O 


ct> 


CM 


^bnincmv  jaiivmjh  viva  j.jj 


M 
I — 
or 

Li) 

X 

o 


o 

CO 


o 

If) 


o 

3- 


o 

ro 


O 

CM 


> 

O 

z 

lx) 

3 

d 

UJ 

X 

u. 


a) 

-t-> 

3 


sz 

-p 
o 
o 
e 
c n 

u 

o 

ft 

>J 

o 

c 

a) 

3 

cr 

0) 

ft 


m 

> 

3 

T3 

3 

■P 


ft 

e 

< 

3 

> 

•rH 

-P 

3 

pH 

3 

cc 

3 

■P 

3 

Q 

Eh 

ft 

ft 


C'- 

cn 

3 

u 

S) 

•H 


68 


27  SEP  78 


00 

r- 

CL 

LU 

<S) 

03 

CM 


— 

CM 

rO 

cf 

m 

X 

X 

X 

X 

X 

0 

O 

0 

O 

CD 

ZD 

ZD 

3 

ZD 

ZD 

O 

0 

O 

O 

O 

X 

X 

X 

CC 

CC 

a 

0 

♦ 

0 

x 

«y 


<M 


**> 


«pO 


<ot«  o 

0*0  o 


(SI 
h- 
cc 
o £ 


O 

CO 


o 

in 


o 

rO 


jO 

CM 


> 

o 

2 

UJ 

ZD 

a 

UJ 

cc 


^ o <o.  cm  oo  cr  o cq 
rO‘  r'i  cm'  c\i  — — ' 

sannidAiv  3 a i iv  13d  viva  133 


CM 


Ih 

O 

U. 

>» 

o 

C 

0 

3 

o* 

3 

(h 

u 


m 

> 

3 

T3  03 
3 .* 
■P  O 
•H  3 
pH  p< 

a a > 
S cc 
< 

o 

<13  03 
> O 
•h  m 
■p 

3 iH 
pH  3 
3 O 
(X  -H 
(h 
3 3 
■P  x: 
3 P. 
a 3 

•H 
6H  g 

Pt.  3 
Ph  X 


ON 

tn 

3 

u 


5> 

•H 


70 


• 12  OCT  78 

* 13  OCT  78 
PLATE2  10  CM 


t 


71 


Figure  40.  Ratio  Of  BWA  To  VWA  vs.  Frequency  (10  cm. 
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Figure  44.  Ratio  Of  BWA  To  VWA  vs.  Frequency  (50  cm. 
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Figure  51  • Ratio  Of  BWA  To  VWA  vs.  Frequency  (40  cm. 
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Figure  60.  Ratio  Of  BWA  To  VWA  vs.  Range  (2?. 5 kHz) 
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Figure  62.  Ratio  Of  BWA  To  VWA  vs.  Frequency  (20  cm. 
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Figure  63.  Ratio  of  BWA  To  VWA  vs.  Frequency  (30  cm. 
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Figure  69.  Ratio  of  BWA  To  VWA  vs.  Frequency  (50  cm. 
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